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In this paper, the lepton flavor violating r~ — >■ fi~ PP{PP — K'^ K~ , K'^ R'^ ,tv'^-k~ ,tv'^tt'^) decays 
are studied in the framework of the two Higgs doublet model(2HDM) III. We calculate these decays 
branching ratios and get the bounds of model parameter jAr^| from the experimental upper limits. 
Our results show that, the neutral Higgs bosons have tree-level contributions to these decays. Among 
these decays, the — > K'^K~ decay is most sensitive to \Xrfi\. In the existing parameters space, 
these decays could reach the measure capability of B factory. These processes can provide some 
valuable information to future research and furthermore present the reliable evidence to test the 
2HDM HI model. 

PACS numbers: 13.35.Dx, 12.15.Mm, 12.60.-i 



I. INTRODUCTION 



Flavor physics have made rapid development in these decades. In addition to B physics, t physics, including 
determination of as from the inclusive hadronic width, charged-current universality tests, and lepton flavor violation 
decays etc., also belongs these days to one branch of flavor physics. Among these topics, lepton flavor violation(LFV) 
decays raises in importance after the discovery of neutrino flavor oscillations and related non-zero neutrino masses [1]. 
In the SM these processes are forbidden or suppressed strongly, therefore, LFV decays could be a sharp tool 
to seek for some new scenarios virith new LFV source and/or new particles. The theoretical investigations of 
T — !• 3Z,T — > /17, T IP{V'^) have sprung up in different contexts [5H5]. With only two meson in final states, 
— > 11^ PP decays are so clean as to provide some information of QCD. And the experimental upper limit of 
T ^ K+K-{K°K°,TT+n-,T:°TT°) are ^6]: 

B^T- fi-K+K-) < 6.8 X 10^^ 90%CL 
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B{t- fi-jT+n-) < 3.3 X 10 



8 



90%CL 



B{t- fi^K'^sK°) < 3.4 X 10-^ 90%Ci 



(1) 



B{t- /i-7r°7r°) < 1.4 x 10 
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90%CL 



(2) 



where the former two values have improved the previous upper bounds by almost an order of magnitude. 

There are also lots of theoretical researches on r IPP decays in many possible extensions of the SMjTHin]. For 
T —5- IPP decays with hadrons in final states, their amplitudes could be separated into leptonic vertexes and hadronic 
parts. One approach to handling the latter is usually parameterized as hadrons mass and their decay constants, 
which could be determined by experimental values, some authors have made analysis of these processes from views 
of vector meson dominance, chiral symmetry breaking, Breit-Wigner propagators, etc.[7j. For the hadronisation of 
final state in t decays, its scale is at the order of IGeV which lies in the non-perturbative region. Hence, we need 
consider the non-perturbative methods one of which is the Chiral Perturbative Theory(xPT)[ll . Different from 
one pseudoscalar meson in final state, the resonances are participated in the processes of — > PP. Stemmed 
from xPT, the Resonance Chiral Theory (.Rx^^) has developed [H]. Using RxT, E. Arganda et al. have investigated 
these processes in two constrained MSSM-seesaw scenarios[8]. M.Herrero et al. also have make an discussion on the 
sensitivity of LFV tau decays the Higgs sector of SUSY-seesaw models[5]. The new particles effects to — > /.i^PP 
decays in the TC model and the LHT model are calculated by Yue chongxing's group [lO]. 

In 2HDM model III, it naturally introduces flavor-changing neutral currents (FCNCs) at tree level. In order 
to satisfy the current experiment constrains, the tree- level FCNCs are suppressed in low-energy experiments for 
the first two generation fermions. While processes concerning with the third generation fermions would be larger. 
These FCNCs with neutral Higgs bosons mediated may produce sizable effects to the t — fi transition. The authors 
in [13] have discussed the t — ^ 3/i decay and the Higgs sector's contributions, r — )■ fiP{V'^) decays have been 
studied under this scenario in our previous workp] where the hadronisation in final state is merely expressed in 
terms of the meson decay constants and meson masses. The r — > fj,P decays could yield one pseudoscalar meson 
from the vacuum state through the scalar and pseudoscalar currents. Hence, this type decay could occur at the 
tree level through the neutral Higgs bosons exchange. In this paper, we extend our discussion to the case of two 
pseudoscalar mesons in the hadronic final state and deal with these four decays by RxT. Our results suggest that, 
the neutral Higgs bosons contribute at the tree level in 2HDM model HI. The model parameter A^/^ is restrained 
at 0(10 ~ lO'^) and the decay branching ratios could as large as the current upper limits. 

The paper is organized as follows: In section II, we make a brief introduction of the theoretical framework for 
the two-Higgs-doublet model HI. In section 111, we briefly introduce the Resonance Chiral Theory. In the next 
section, we deliberate the calculation of the decay amplitudes with Resonance Chiral Theory and our numerical 



predictions. Our conclusions are listed in the last section. 



II. THE TWO-HIGGS -DOUBLET MODEL III 

As the simplest extension of the SM, the Two-Higgs-Doublet Model has an additional Higgs doublet. In order 
to ensure the forbidden FCNCs at tree level, it requires either the same doublet couple to the w-type and rf-type 
quarks(2HDM I) or one scalar doublet couple to the w-type quarks and the other to d-type quarks(2HDM II). 
While in the 2HDM III[T4l[T5], two Higgs doublets could couple to the w-type and rf-type quarks simultaneously. 
Particularly, without an ad hoc discrete symmetry exerted, this model permits flavor changing neutral currents 
occur at the tree level. 

The Yukawa Lagrangian is generally expressed as the following form: 



(3) 



where Hi{i = 1,2) are the two Higgs doublets. Qi^L is the left-handed fermion doublet, Uj^r and Dj ji are the 
right-handed singlets, respectively. These Qi,L, Uj^R and Dj R are weak eigenstates, which can be rotated into mass 
eigenstates. While rj^'^ and ^'^■^ are the non-diagonal matrices of the Yukawa couplings. 
We can conveniently choose a suitable basis to denote Hi and H2 as: 

/ 



Hi 
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71 




(4) 
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where G°'^ are the Goldstone bosons, and are the physical charged-Higgs boson and CP-odd neutral Higgs 
boson, respectively. Its virtue is the first doublet Hi corresponds to the scalar doublet of the SM while the new 
Higgs fields arise from the second doublet H2- 

The CP-even neutral Higgs boson mass eigenstates iJ" and /i" are linear combinations of (pl and 4>2 i'^ Eq.Q, 



H° = 05 cos a -I- (/)2 sin a, ft," = sin a -I- 0^ 



cos a, 



(5) 



where a is the mixing angle. 

After diagonalizing the mass matrix of the fermion fields, the Yukawa Lagrangian becomes|16j 
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X~DV^KM [MuR - MdL] U ~ —x^UVcKM [MdR ~ MuL\ D 



\u i^R-i^^L U + D i^^L-i^R 



—;=U \ — Af[/cosa + 
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sina > [/ < — Af^ cos a + 
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(6) 



where U and D now are the fermion mass eigenstates and 



U.D-^-l JJ^D yU,D 



V — V L J 'I "R. 



(7) 
(8) 



where V^'^^ are the rotation matrices acting on up and down-type quarks, with left and right chirahties respectively. 
Thus VcKM ~ i^Ey^F is usual Cabibbo-Kobayashi-Maskawa (CKM) matrix. In general, the matrices i)^'^ 
of Eq.l?]) are diagonal, while the matrices ^^'^ are non-diagonal which could induce scalar-mediated FCNC. Seen 
from Eq.Q, the coupling of neutral Higgs bosons to the fermions could generate FCNC parts. For the arbitrariness 
of definition for ■fj^'^ couplings, we can adopt the rotated couplings expressed ^^'^ in stead of hereafter. 



In this work, we use the Cheng-Sher ansatz[l4] 



U.D _ . V2^/m-r 



V 



(9) 



which ensures that the FCNCs within the first two generations are naturally suppressed by small fermions masses. 
This ansatz suggests that LFV couplings involving the electron are suppressed, while LFV transitions involving 
muon and tau are much less suppressed and may lead to some loop effects which are promising to be tested by the 
future B factory experiments. In Eq.([9]), the parameter Xij is complex and i,j are the generation indexes. In this 
study, we shall discuss the phenomenological applications of the type III 2IIDM. 



III. THE RESONANCE CHIRAL THEORY 

For the intermediate and low energy parts of hadronic spectrum, they locate at the non-perturbative region 
where ordinary perturbative QCD methods does not work. Hence, many non-perturbative approaches have been 
developed, such as xPT llj, QCD sum rules [1 7 , lattice gauge theory[18] and so on. The large-A'c expansion of 
SU{Nc) QCD[Tn] is a suitable idea. 

xPT is a appropriate method with 1/Nc expansion which is very successful in the energy reign of ~ IGeV. It 
could deal with — > /i^P decays which will be discussed in our later paper. While for — >■ /i^PP processes, the 
resonances paly a dynamical role so that we should by the aid of RxT^2\- Motivated by partially and large 
Nc QCD, RxT could describe the immediate energy region. Its advantages are not only to realize the non-linear 



of spontaneous chiral symmetry breaking, but also to study meson structure without imposing any structure in 
advance. Some application on i?xT have been elaborated in [51121 HOI- For energies 1 ^ 2GeV, we restrict it to 
only the lightest resonance in each channel. The detailed deliberation can be found in Ref. ^ . 

The symmetry of QCD breaks spontaneously from SU{3)R^ SU{3)l to SU{3)v and produces eight Goldstone 
bosons in the spectrum. We regard these Goldstone bosons as the lightest hadrons. The xPT Lagrangian can be 
constructed from pseudoscalar fields and external source w*j(a;), a^a;), Si(a;),p*(a;): 



C — £o + £i, 

= Y' « =<y' S^SiXi, p^piXi^ 



(10) 



where the Ci is the massless QCD Lagrangian. G^^ denotes the gluon fields, v^, a'^, s,p are matrices in the flavor 
fields and A^s are the Gell-Mann matrices. The QCD generating functional Zqcd[v, a, s,p] could be written: 



JZQco[v,a,s,p] 



[DG,,][Dq][Dq]e 



(11) 



We introduce the lightest U{3) nonet of pseudoscalar mesons: 



^^^^ = 



a=0 



( 1 ^0 I 1 
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The unitary 3x3 matrix u{x) can be presented as: 



76' 
^0 



76^ ^ 73^ / 



(12) 



(13) 



The leading 0{j?-)xPTSU['i)L x S'C/(3)fl chiral Lagrangian is 



F 



= i[u^((9p - irp)u - u(5p - iZ^)-u^], x+ = u^x"^ + ^^X^", X = 2Bo(s + ip). 



(14) 



Interactions with electroweak bosons can be accommodated through the vector = (t^i+'h) 



and axial-vector 



external fields. The effective coupling constant F is approximately equal the pion decay constant of 



TT. In the chiral limit there has BqF^ = — (O|'i/;'0|O)o. The chiral tensor i/j provides masses to the Goldstone bosons 



/ 



through the external scalar field s ^ Ai + ■ ■ ■ ,A4 
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Indeed in the isospin limit we have: 
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Boniu = Bound = 7:m^, Bonis = m 
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(15) 



The mass eigenstates 77 and 77' are related to the octet T^g and singlet 770 states through the rotation matrix: 



COS f — sm f 



sm V cos ( 



1 









(16) 



where 6* = 18". 

In the processes of r decaying to two pseudo-scalars, vector and scalar resonances generally play a propelling 
role. However, due to the higher masses of scalar resonances, their effects are ignored, 8 . Then we carry out 
antisymmetric tensor fields to introduce vector resonances. The nonet of resonance fields V^i/ reads: 
/ 1 „o _i_ 1 1 \ 



-T2P +71^8 + 73^0 



73 



tJ3 



V 



K* 



K 
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(17) 



Then the i?xT Lagrangian puts: 



RxT 



+ '^f2'|: - £fcm + 

1,_. M? 



(2)' -^V/ — '-km ^ -^(2)' 
/1//-2 



(18) 



where Cv is the resonance Lagrangian. The relevant definitions can be found in Ref.[S|. The corresponding QCD 
generating functional is written as 



piZQCD[v,a,s,p\ 



[Du][DV]e 



i J d'^xCn.xTi'f'jVjV^a^Sjp] 



(19) 



Through making the proper partial derivatives of the functional action, we can get the hadronisation of bilinear 
quark currents: 



The final results from Eq. ( 20 ) are 



f"^ f 



(20) 



(21) 



/I* 



\boF^ < X'{u^u^ +uu) > 
^BqF'^ < A'(m^u^ - uu) > 



(22) 
(23) 
(24) 



IV. THE DISCUSSION FOR r ^ ^ PP DECAYS 



In 2HDM model III, the neutral Higgs bosons mediated tree diagrams have contributions to fi^PP 
processes. The amplitudes could be factorized into leptonic vertex corrections and hadronic parts described with 
hadronic matrix elements, which express as: 

iGp 



{fl-PP\M\T- 



2V2 



H"* ■ Xr^ ■ iMs+p + H<i ■ X* ■ {fir)s-p 



<PP\iqq)s\0> 



N''* ■ ■ {fir)s+p - N'^ ■ ■ (M 



S-P 



<PP\{qq)p\0> 



(25) 



where H'^^N'^ are auxiliary functions can be find in Appendix. From the Eq.(25l, we can see that three neutral 



Higgs bosons perform roles to two pseudo-scalar mesons through {qq)s±p operators. It differs from the case of 
MSSM models where the 7 contributions is the dominate one and only i/" and take effects at the large tan /3[8]. 
For the heavy Higgs bosons, the hadronic final state is not sensitive to resonances and known little. It should be 
noted that the pseudo-scalars currents have contributions to only one pseudo-scalar meson in final states. So we 



have dealt with the hadronic matrix elements by virtue of Eq.(|23|),(|24|) and the following currents: 

1 

1 

2' 
I 



uu = -S-" 



2^3 03 



-dd 



2 2^3 V6 



-ss 



V3 Ve 



The obtained amplitudes read as: 



M{t- H^PP) = 



iGp 
2a/2 



y^TV^ • [T{PP) ■ Xr^ ■ {flT)s+P + T*{PP) ■ X*^ ■ iMs-p] 



(26) 



(27) 



where 



T{K+K-) = A ■ [Re{Xuu) ■ ml + i?e(A,,) • {2m\ - ml)] + Bi ■ [Im{Xss) ■ (2m^ - ml) - /m(A„„) • ml] 

+F ■ m% (28) 
TiK°K°) = A ■ [Re{Xdd) ■ ml + Re{Xss) ■ {2m\ ~ ml)] + Bi ■ [/m(A,,) • (2m^ - ml) + ImiXdd) ■ ml] 

+F ■ m% (29) 
T{tt+tt-) = ml ■ {A[Re{Xuu) + Re{Xdd)] + Bi ■ [Im{Xdd) ~ /m(A„„)] + F} (30) 



T(7rV) 
A 



2^/2 

. 2 

sm a 



{A[i?e(A, 



''HO 



cos^ a 



1 



(31) 



F — 2sinQ;cosQ!(- 



•-HO 



TABLE L Constraints on the \ij in quark and lepton sector. 



Quark sector 


Bounds and restrictions 


Process and Restriction 


References 


|A„„|,|Addl ~ 0(1) 


F» - FO mixing (F ^K,Bd, D),Rb, P,B ^ X.-y 


El 


|A„| = 0.02, |Af,5| = 50 


Bd — Bd mixing ,6 — > 37 


M 


|A„| < 0.5, |A6i,| < 70, \\tt\bb\ ~ 3, |A,,1 G [80, 120] 


Bd - Bd mixing ,b ^ s-y, po, Rb, NEDM 


m 


Att = Ate = 0.1, \\bb\ = Afts = 50 


h" ^ ,// 


m 


\\tt\ — 0.3, \Xbb\ = 35, Xij = 


B PP, PV, VV 


m 




Arr = Ap^ = 5,50, 


Bd.s i+r 


m 


A^i^ = Arr = A^r = Ae^t = 10 


^ ,// 


m 


A,^ ~ 0(10) -0(10") 


{g — 2)^,m_40 — > 00 


m 


A,^~O(102)-O(10^) 


r — > 3/i, r — > /i7 


m 



^ Note that the constraints in[26] and [27] are denoted by our notation. 



V. NUMERICAL RESULTS 



In our calculation, the input parameters are the Higgs masses, mixing angle a, and their phase angles 9ij. 
We using the values of neutral Higgs boson masses in literature [22l [231 [26], where the experimental constraints of 
B — B mixing, b — > 57, p", Ri, considered. 

rriHO = leOGeV, m^o = llbGeV, m^o = UOGeV, a = 7r/4, 6 = n/4 (32) 

According the mesons quark constants in final states, the involved factors of quark sector are A„„, Xdd and Ass. The 
bounds of |A,-^| from different phenomenological considerations [231 [13 [IZ] are demonstrated in Tab. I, too. For the 
first generation FC couplings are suppressed, the values of Xuu and Xm are less than 1 21J. And the Bq — Bq mixing 
constrains approximately Ass in 80 ^ 120|23|. These bounds are considered in our calculation. In the following 
paragraphs, we will analysis the relations of these decays branching ratios and relevant parameters . First, we take 
lAr/il = 5 and study the relations of branching ratio and other parameters. 

The Fig.l shows Br for r~ /.i^KK decays versus model parameters, where the left figure is Br{T^ — 
fi^ K^K^) versus |A,j„| and |Ass|, and the right is Br{T^ ijTK'^K^) versus \Xdd\ and |Ass|. One can see 
from Fig. 1(a) that Br{T^ — > fi^ K^) raises with the increase of |Ass| but does not vary with |A„„| growing. 
For ^^K^K^ decay, the same as that of ji^K^K^ decay except that \Xdd\ replaces |A„„|- In the 



1. 2 

Br . 




FIG. 1: Left:i3r for r — )■ M K'^ K versus \\uu\ and lAss|; Right:_Br for r — >■ /i K^K'^ versus \\dd\ and |Ass| with 
|A,^| = 5. 




FIG. 3: Left: Br for r~ — >■ ^~KK versus |At,j| for |Auu| — 0.5, |Ass| = 100. The sohd hue stands for r~ — >■ K~ , the 

dashing hne for — > K"K^, the horizon hne for current experimental upper limit for — >■ fi^ K~ . Right:_Br for 

— > fi^TVTT versus |At-^| for |A„u| = \\dd\ = 0.5. The solid hne stands for — >■ /i^Tr^vr", the dashing line for — >■ /i~7r''7r'^, 
the horizon line for current experimental upper limit for — >■ h^tt^tv^ . 



mentioned parameter spaces, both two decays ratios could reach the order of 10^ — 8). Comparing Eq.(24) and 
(25), we could find that the structures of these decay amplitudes are similar. Hence Br{T~ jjTK^K^) also rises 
with the increase of \\as\ but is almost not affected by the modification of \\dd\- It resulted from the suppressed 
\^uu\^ \^dd\ and that \\ss\ is larger than |A„„|, \ \dd\ two order of magnitudes. 

The functions of r~ ^~'inr decays versus model parameters are presented in Fig.2, where the left figure is 
Br{T~ /Lt~7r+7r~) versus |A„„| and |Arfd|, and the right is Br{T~ /i~7r°7r'^) versus |A„„| and \Xdd\- Both 
decay amplitudes are relevant to |A„„| and \Xdd\- When |A^^| = 5, their decay ratios move upward with |At,„| 
and \Xdd\ and both extend the order of 10^ — 14). We find that Br{T~ /i~7r+7r~) climbs more rapidly than 
Br{T~ /i~7r°7r'^). It is because the contributions of |A„„| and \Xdd\ to the former amplitude are larger than those 
to the latter amplitude. 

Then, we take |A„„| = \Xdd\ = 0.5, |Ass| = 100 and analysis the relation of branching ratios versus \Xrij.\. Fig. 3 
gives four decays branching ratios versus |At-^| for other fixed parameters. The left figure is branching ratios 
for T~ — )• ^~KK, where the solid line stands for r~ — )• n~K~^K~, the dashing line for r~ — )• fjTK^K^. The 
experimental data for r~ — >■ iJ,~K~^K~ and t~ — > ij,~K°K° are denoted by the lower horizon line and upper 
horizon dash line, respectively. From left figure, we could see the lower horizon line constraints |At^| at ~ 4.5 
for Br{T~ — )• ii~ K^K~), and the upper limit for Br{T~ fjTK^K^) constraints |At-^| at the order of 0(1). It 
should be noted that the latest experimental value of Br{T~ — )• ^~K~^K~) is based on 671 fb~^ data. If the 
data of Br{T~ — >■ n~K^K^) have been updated, the bound obtained will be more respective. The right figure is 
Br for T~ iJ.~7r'jT where the solid line stands for r~ — )• /i~7r"'"7r~ and the dashing line for r~ /x~7r°7r°. The 
experimental data for r~ fj.~iT~^n~ and t~ /U~7r°7r'^ are denoted by the lower horizon line and upper horizon 
dash line, respectively. From the right figure, we could find Br(r~ — >■ /i~7r+7r~) is more sensitive to |At-^| than 
Br{T~ — >■ /Lt~7r°7r°). In all, the r~ — >■ ijl~K^K~ decay would make tighter constraints on the Higgs couplings than 
those from decays. 



VI. CONCLUSION 

Sum up, we have calculated the branching ratios of t~ — >■ ^,~PP{PP = K~^K~ , K^K'^ ,n'^Tr~ ,tt^tt°) decays in 
the model III 2HDM. The neutral Higgs bosons contribute to theses decays at the tree-level. Comparing to the 
T~ fj.~P decays, the resonances play a part in t~ — )• n~PP processes and to which the massive Higgs bosons 
have insensitivity. Only the scalar currents contribute to these decays. Our work suggests that the parameter |Ax|u| 
is constrained at the order of 0(1 ~ 10^) by the experimental data. And in the rational parameters space, their 
branching ratios can reach the experimental values. The r~ — ^ n~K~^K~ decay is most sensitive to |At^|. Our 



study is hoped to supply good information for the future experiment and explore the structure of the 2HDA/[ III 
model. 
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